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THE METRIC-RESTRICTED INVERSE DESIGN PROBLEM 


AMIT ACHARYA, MARTA LEWICKA AND MOHAMMAD REZA PAKZAD 

Abstract. We study a class of design problems in solid mechanics, leading to a variation on the 
classical question of equi-dimensional embeddability of Riemannian manifolds. In this general new 
context, we derive a necessary and sufficient existence condition, given through a system of total 
differential equations, and discuss its integrability. In the classical context, the same approach 
yields conditions of immersibility of a given metric in terms of the Riemann curvature tensor. 
In the present situation, the equations do not close in a straightforward manner, and successive 
differentiation of the compatibility conditions leads to a new algebraic description of integrability. 
We also recast the problem in a variational setting and analyze the infimum of the appropriate 
incompatibility energy, resembling the non-Euclidean elasticity. We then derive a T-convergence 
result for dimension reduction from 3d to 2d in the Kirchlioff energy scaling regime. 


1. The metric-restricted inverse design problem 

Assume T is a manifold of material types, differentiated by their structure, density, swelling- 
shrinkage rates and other qualities. We let any material type T £ T be naturally endowed with a 
prestrain g(T), where g : T —>• Ksym, P os is a given smooth mapping, taking values in the symmetric 
positive definite tensors. 

Suppose now that we need to manufacture a 2-dimensional membrane S C K 3 , where at any 
given point p £ S a material of type T(p) £ T must be used for a given T : S —> T- The question 
is how to print a thin film [/Cl 2 in a manner that the activation u : U —> M 3 of the prestrain in 
the film would result in a deformation leading eventually to the desired surface shape S. 

The above described problem is natural as a design question in various areas of solid mechanics, 
even though the involved tensors are not intrinsic geometric objects. For example, it includes the 
subproblems and extensions to higher dimensions: 

(a) Given the deformed configuration of an elastic 2-dimensional membrane and the rect¬ 
angular Cartesian components of the Right Cauchy-Green tensor field of a deformation, 
mapping a flat undeformed reference of the membrane to it, find the flat reference config¬ 
uration and the deformation of the membrarifj]] 

(b) Given the deformed configuration of a 3-dimensional body and the rectangular Carte¬ 
sian components of the Right Cauchy-Green tensor field of the deformation, mapping a 
reference configuration to it, find the reference configuration and the deformation. 

(c) Suppose the current configuration of a 3-dimensional, plastically deformed body is given, 
and on it is specified the rectangular Cartesian components of a plastic distortion F p . Find 
a reference configuration and a deformation £, mapping this reference to the given current 
configuration, such that the latter is stress-free. Assume that the stress response of the 
material is such that the stress vanishes if and only if ('V((F p )~ * 1 ) T (V ((Fp)^ 1 ) = M 3 . 


Date : today. 

■'■We thank Kaushik Bhattacharya for bringing this problem to our attention. 
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In view of natal [19], the activation u must be an isometric immersion of the Riemannian 
manifold ( U,Q ) into R 3 , where Q is the prestrain in the flat (referential) thin film. Our design 
problem requires hence that we find an unknown reference configuration U C R 2 , an unknown 
material distribution T : U —> T and an unknown deformation u : U R 3 such that 

(i) S = u(U), 

(ii) For any x E U, the point u(x) carries a material of type T(u(x)), i.e. T(u(x)) = T(x), 

(iii) 'Vu(x) T 'Vu(x) = G(x) := g(fZ(x)). 

If the membrane 5 C R 3 is a smooth surface, then letting g := goT : S —> M 2 ^ pos , the conditions 
(i)-(iii) simplify to finding a domain U Cl 2 and a bijection u : U S, such that: 

(1.1) (Vrx) T (Vu)(x) = g(u(x)) Vx 6 U. 

The smoothness of g is determined by the regularity of T and of the mappings g and T. 


Essentially, in all of the applications defined above (e.g. membrane, 3-d), we are dealing with a 
general class of nonlinear elastic constitutive assumptions involving pre-strain, with the require¬ 
ment that the stored energy density evaluated at the Identity tensor (of appropriate dimension¬ 
ality) attain the value zero; in other words, we look for stress-free deformations of a prestrained 
body. Given a prestrain field specified on the target configuration , we explore the question of exis¬ 
tence of deformations that allow such a minimum energy state to be attained pointwise, as well as 
the characterization of the constraints on the pre-strain field that allows such attainment. In the 
language of mechanics, note that the question ( 1 . 1 ) may be rephrased as looking for deformations 
u of the reference U such that: 


Vu 


--1 


Vr 


u) 


-i 


= Id, 


where the expression on the left-hand-side of the above equality is the sole argument of the frame- 
indifferent nonlinear elastic energy density function of the material. Thus, \Jg(u) needs to be 
capable of annihilation by the right stretch tensor of a deformation, a condition that is expressed 
in terms of spatial derivatives of g{u) on U; the main difficulty is that both u and U are unknown, 
so this differential condition cannot simply be written down and a more sophisticated idea than 
the standard vanishing of the Riemann-Christoffel tensor of a metric is needed. 


The inverse design problem that we study in this paper, can be further rephrased as follows. 
Let y : 17 —» R 3 , be a smooth parametrization of S' = y(fi). Find a change of variable £ : —> U 

so that the pull back of the Euclidean metric on S through y is realized by the following formula: 

( 1 . 2 ) (Vy) T Vy = (V £) T (5 ° y)V£ in IL 


Clearly, once a solution £ of ( 1 . 2 ) is found, the material type distribution T, which is needed for 
the construction of the printed film U, can be calculated by: 

■= T o y o £ -1 : U = £(f2) —>• T, 


since u = y o £ 1 satisfies (1.1) and consequently the properties (i)-(iii) hold. Any £ satisfying 

(1.2) is an isometry between the Riemannian manifolds and (U,G o £~ x ), with metrics: 

(1.3) G = (X7y) T Vy and G = goy on 


n. 


For convenience of the reader, we gather some of our notational symbols in Figure 0 

The same problem can be set up for a three dimensional shape S C R 3 . In that case, the pre¬ 
strain mapping g must take values in R sym,pos an ^ the printed prestrained reference configuration 
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is modeled by U C M . The equation to be solved is still given by (1.2), now posed in C M 3 . 
Equivalently, a solution u : U S to (1.1) is obtained as in (iii) above for T := T o u and it is 
the absolute minimizer of the prestrain elastic energy (see e.g. 


(1.4) E{u, U, Q) := / dist 2 ((Vu)^ -1 / 2 , 50(3)) dx = 0. 

Ju 

Here, dist(F, 50(3)) stands for the distance of a matrix F from the compact set 50(3), with 
respect to the Hilbert-Schmidt norm. What distinguishes our problems from the classical isometric 
immersion problem in differential geometry, where one looks for an isometric mapping between 
two given manifolds (fi,G) and (U,G), is that the target manifold U = £(fl) and its Riemannian 
metric Q = Go( _1 are only given a-posteriori, after the solution is found. Note that only when G 
is constant, the target metric becomes a-priori well defined and can be extended over the whole 
of M n , as it is independent of £, and then the problem reduces to the classical case (see Example 



Figure 1.1. Geometry of the problem. 


Remark 1.1. Note that the minimizing solution in the equidimensional problem (1.4) is unique 
up to rigid motions. Hence, If the configuration (U, Q) is printed, the elastic body is bound to take 
the required shape 5 as requested in the design problem. Note that in case the required shape 5 
is two dimensional, uniqueness fails due to one more degree of freedom as the deformation u of 
[/Cl 2 takes values in a higher dimensional space M 3 . In this case, other restrictions on u have 
to be imposed, to solve the original design problem. One particular remedy is to consider the 
membrane as a thin three dimensional body; we implement this approach in Section [3j Note that 
solving (1.1) directly in two dimensions implies that a compatible reference configuration can be 
found, as formulated in (a) above. 

Another approach in the two dimensional case would be to use a parametrization y : H —> S 
which is bending energy minimizing among all other mappings y : fl —> 5 that induce the same 
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metric (Vy) r Vy. This would necessitate a study of the bending energy effects, which is beyond 
the scope of the present paper. For the sake of comparison, we mention a parallel problem, studied 
in |15| , where the authors calculate the prestrain to be printed in the thin film by minimizing a full 
three dimensional energy of deformation, consisting of both stretching and bending. The nature 
of the problem in m is different than ours, in as much as the material constraints are not present. 


The paper is organized as follows. In Sections 2 and 3, we study the above mentioned variational 


formulation of (1.2) and analyze the infimum value of the appropriate incompatibility energy, 


resembling the non-Euclidean elasticity m • We derive a T-convergence result for the dimension 
reduction from 3d to 2d, in the Kirchhoff-like energy scaling regime, corresponding to the square 


of thickness of the thin film. In Sections 4-9, we formulate (1.2) as an algebraically constrained 


system of total differential equations, in which the second derivatives of £ are expressed in terms 
of its first derivatives and the Christoffel symbols of the involved metrics. The idea is then 
to investigate the integrability conditions of this system. When this method is applied in the 
context of the standard Riemannian isometric immersion problem, the parameters involving £ 
can be removed from the conditions and the intrinsic conditions of immersibility will be given 
in terms of the Riemann curvature tensors. In our situation, the equations do not close in a 
straightforward manner, and successive differentiation of the compatibility conditions leads to a 
more sophisticated algebraic description of solvability. This approach has been adapted in [1] for 
deriving compatibility conditions for the Left Cauchy-Green tensor. 

Acknowledgments. This project is based upon work supported by, among others, the Na¬ 
tional Science Foundation. A.A. acknowledges support in part from grants NSF-CMMI-1435624, 
NSF-DMS-1434734, and ARO W911NF-15-1-0239. M.L. was partially supported by the NSF 
grants DMS-0846996 and DMS-1406730. M.R.P. was partially supported by the NSF grant DMS- 
1210258. A part of this work was completed while the second and the third authors visited the 
Forschungsinstitut fur Mathematik at ETH (Zurich, Switzerland). The institute’s hospitality is 
gratefully acknowledged. 


2. A VARIATIONAL REFORMULATION OF THE PROBLEM (1.3) 


In this section, we recast the problem (1.3) in a variational setting, similar to that of non- 


Euclidean elasticity (19] . Using the same arguments as in mm, we will analyze the properties of 
the infimum value of the appropriate incompatibility energy, over the natural class of deformations 
of IF 1 ' 2 regularity. 


We will first discuss the problem (1.2) in the general n-dimensional setting and only later restrict 
to the case n = 2 (or n = 3). Hence, we assume that fl is an open, bounded, simply connected 
and smooth subset of M n . We look for a bilipschitz map £ : —> U := £(fi), satisfying (1.3) and 

which is orientation preserving: 

(2.1) detV£>0 in fL 


We begin by rewriting (1.3) as: 

( 2 . 2 ) 


G = (G 1 / 2 VC) T (G 1 / 2 VC). 

Note that, in view of the polar decomposition theorem of matrices, a vector field £ : H 


is a 


solution to (2.2), augmented by the constraint (2.1), both valid a.e. in Q, if and only if: 
(2.3) Va.e. iGll 3R = R(x) G SO{n ) G 1/2 V£ = R G 1/2 , 
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where G 102 (x) denotes the unique symmetric positive definite square root of G(x) E 
while SO(n) stands for the set of special orthogonal matrices. Define: 


»nxn 
^ sym,pos ’ 


(2.4) 


F(£) = / dist 2 (G 1 / 2 (V^)G _1/2 ,S'0(n)) dx V£ E W / ^’ c 1 (n,M n ), 

Jn 


where dist (F,SO(n)) is the calculated distance of a matrix F from the compact set SO(n), with 
respect to the Hilbert-Schmidt norm of matrices. It immediately follows that F(£) = 0 if and only 
if £ is a solution to \2.2\ and hence to (1.3), together with (2.1). Also, note that E {£) < oo if and 
only if £ E VF 1 ’ 2 (fl,M n ), as can be easily deduced from the inequality: 


(2.5) 


VF E 


|F| 2 < C\G 1/2 FG~ 1/2 \ 


valid with a constant C > 0 independent of x and F. 

Finally, observe that, due to the uniform positive definiteness of the matrix field G: 


( 2 . 6 ) 


Iff = £ |Fej| 2 < C^ j (Fe i ,GFe i ) < Ctrace ( F T GF ). 


i= 1 


i=1 


Proposition 2.1. (i) Assume that the metrics G,G are C(Ll,W nxn ) regular. Let £ E W^ 1 (D, R n ) 
satisfy (0 for a.e. Then £ E W 1,0O (fi, R n ) must be Lipschitz continuous. 

_ (ii) Assume additionally that for some k > 0 and 0 < /i < 1, G,G E C fc,/i (D, M nXTl ). If (1.3), 

(2.1) hold a.e. in Q (so that F(£) = 0), then £ E C fc+ 1 ,/i (D,M n ). 

Proof. The first assertion clearly follows from the boundedness of G and positive definiteness of 
G, through ( |2.6[ ). To prove (ii), recall that for a matrix F E M nxn , the matrix of cofactors of F is 
cof F, with (cof F)ij = (—l) i+J det Fjj, where Fij E is obtained from F by deleting 

its ith row and jth column. Then, (1.3) implies that: 

detV£ = (^^) 1/2 = : a E C(D,M+) and cofV£ = aG(V£)G -1 . 

Since div(cof V£) = 0 for £ E W 1,oa (where the divergence of the cofactor matrix is always taken 
row-wise), we obtain that £ satisfies the following linear system of differential equations, in the 
weak sense: 

div(aG(V£)G" 1 ) = 0. 

Writing in coordinates £ = (£ x .. •£ ri ), and using the Einstein summation convention, the above 
system reads: 

Vi = 1 ... n d a (aGijG 0 a d^ j ) = 0. 

The regularity result is now an immediate consequence of m Theorem 3.3] in view of the ellipticity 
of the coefficient matrix Af 0 = aGijG a0 . ■ 

We now prove two further auxiliary results. 

Lemma 2.2. There exist constants C,M > 0, depending only on ||G||l°o and ||G||/,oo ; such that 
for every £ E lF 1 , 2 (D,M n ) there exists £ E lF 1 , 2 (D,M n ) with the properties: 

||V£||l°o < M, ||V£ - veill 2 (n) <GF(£) and F(£) < GF(£). 
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Proof. Use the approximation result of Proposition A.l. in [HI to obtain the truncation £ = £ A , 
for A > 0 having the property that if a matrix F E M nxn satisfies \F\ > A then: 

|F | 2 < Cdist 2 {G 1/2 FG~ 1/2 (x), SO(n )) Vx E SI. 

Then ||V£ a ||l°° < CX := M and further, since V£ = V£ A a.e. in the set {|V£| < A}: 


l|V?-V{ A ||^ (n) = 


'{|VC|>A} 


|V £| 2 < c 


{I V£|>A} 


dist 2 (G 1 / 2 V£G- 1 / 2 ,SO(n)) dx < G£(£). 


The last inequality of the lemma follows from the above by the triangle inequality. 

Lemma 2.3. Let £ E IU 1 , 00 (fi, M n ). Then there exists a unique weak solution cf : n 

j di v{aG{S/(f)G~ 1 ) =0 in SI, 


to: 


(2.7) 


= £ on dSl. 


Moreover, there is constant C > 0, depending only on G and G, and (in a nondecreasing manner) 
on ||V£|| 

Z/°° j such that: 

l|V(£ - 0)||| 2(n) <C£(£). 

Proof. Consider the functional: 

I(<p):= [ (G(V<p)G~ 1 (x) :Vip(x)) dx = [ |o 1 / 2 G 1 / 2 (V^)G '- 1 / 2 | 2 dx Xhp E fU 1 , 2 (n,M n ). 

Jfi Jfi 

The formula (2.5), in which we have implicitly used the coercivity of G and G, implies that: 

||V^||| 2(n) < Cl {ip). 

Therefore, in view of the strict convexity of /, the direct method of calculus of variations implies 
that I admits a unique critical point cf in the set: 

{p E lU 1 , 2 (n,M n ); :p = £ on 


By the symmetry of G and G , ( |2.7[ ) is precisely the Euler-Lagrange equation of I, and hence it is 
satisfied, in the weak sense, by cf. 


Further, for the correction if = £ — (f E IUQ’ 2 (f7,M n ) it follows that: 


V 77 E Wq’ 2 (SI, M n ) 


aG(S/if)G~ l : V?? dx = / aG(V£)G“ i : Vr? - / aG(V</>)G _i : V?? 

J n J q 

= [ aG(V£)G _1 : Vr? 

Jn 

= [ aG(V£)G _1 : V 77 — f cofV£ : V 77 . 

J J q 

Indeed, the last term above equals to 0, since the row-wise divergence of the cofactor matrix of 
V£ is 0, in view of £ being Lipschitz continuous. Use now 77 = if to obtain: 

||VY>||| 2(n) < Cl {if) = cJ iaGiVQG- 1 - cofV£) : S/if dx 

_! 2\V 2 


< C'||V^|U 2(n ) aG^G " 1 - cofV£ 

< C\\Vif\\ L , {n) E^) 1/2 . 
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The last inequality above follows from: 

V|F| < M Vx&n \aGFG~\x) - cofF | 2 < C M dist 2 (G 1/2 FG~ 1/2 ,SO(n)), 
because when G 1 / 2 FG~ 1 / 2 E SO{n) then the difference in the left hand side above equals 0. 


Theorem 2.4. Assume that the metrics G, G E 1 


are Lipschitz continuous. Define: 


( 2 . 8 ) 


k(G,G)= inf E(£). 


Then, k(G,G ) = 0 if and only if there exists a minimizer £ E W 1,2 (fi,M n ) with E(tf) = 0. In 
particular, in view of Proposition 2.1, this is equivalent to £ being a solution to (1.3) (2.1), and £ 
is smooth if G and G are smooth. 

Proof. Assume, by contradiction, that for some sequence of deformations £& E W 1 , 2 (J1, M”), there 
holds limfc^oo E(£f.) = 0. By Lemma 2.2 replacing £*, by £*,, we may without loss of generality 
request that ||V£fc||L°° < M. 

The uniform boundedness of V£fc implies, via the Poincare inequality, and after a modification 
by a constant and passing to a subsequence, if necessary: 

(2.9) 


lim £ fc = £ 
k—> oo 


weakly in W 1,2 (D) 


Consider the decomposition £*. = 4>k + wher e <f>k solves dTfl ) with the boundary data fik = £fc 
on dfi. By the Poincare inequality, Lemma 2.3 implies for the sequence ipk E TTq’ 2 ^): 


lim ipk = 0 strongly in W 1,2 (Q). 

k—>oo 


In view of the convergence in (2.9), the sequence (f>k must be uniformly bounded in VP 1 , 2 (fl), and 


hence by (23J Theorem 4.11, estimate (4.18)]: 

vn' cc n vfc Wf’kWw 2 ’ 2 ^ 1 ) ^ ^a'll^fclliv 1 - 2 ^) < c. 

Consequently, fik converge to strongly in Wl o 2 (IX). Recalling that E(£k) converge to 0, we finally 
conclude that: 

This proves the claimed result. 


m) = o. 


3. A DIMENSION REDUCTION RESULT FOR THE ENERGIES (2.4) 


The variational problem induced by (2.4) is difficult due to the lack of convexity. One way 


of reducing the complexity of this problem is to assume that the target shapes are thin bodies, 
described by a “thin limit” residual theory, whic h is potentially easier to analyze. We concentrate 
on the case n = 3 and the energy functional (2.4) relative to a family of thin films D h = uj x (— |, |) 
with the midplate w CK 2 given by an open, smooth and bounded set. This set-up corresponds to a 
scenario where a target 3-dimensional thin shell is to be manufactured, rather than a 2-dimensional 
surface as in Figure 0 As we shall see below, an approximate realization of the ideal thin shell 
is delivered by solving the variational problem in the limit of the vanishing thickness h. 

Assume further that the given smooth metrics G, G : Cl h —> M 3x3 are thickness-independent: 


G,G(x',x 3 ) = G,G(x') 


Vx = (x',x 3 ) E u x (-^,^), 
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and denote: 
(3.1) 


1 


E h (t h ) = - I W(G 1/2 (Vt h )G~ 1/2 ) dx Vf h G W ll 2 (n h ,R 3 ). 


in h 


The energy density W : M 3x3 —> M + is assumed to be C 2 regular close to 50(3), and to satisfy 
the conditions of normalisation, frame invariance and bound from below: 

3c > 0 VFG M 3x3 VR G 50(3) W(R) = 0, W(RF) = W(F), 

W(F) > cdist 2 (F,50(3)). 


Following the approach of HU , which has been further developed in H2 CE1 El HHI (see also 
for an extensive review of the literature), we obtain the following T-convergence results, describing 
in a rigorous manner the asymptotic behavior of the approximate minimizers of the energy (3.1). 

Theorem 3.1. For a given sequence of deformations £, h G W l,2 (Q h , M 3 ) satisfying: 

(3.2) 30 > 0 V/i E h (f h ) < Ch 2 , 

there exists a sequence of vectors c h G M 3 , such that the following properties hold for the normalised 
deformations y h G lF 1,2 (fl 1 ,M 3 ): 

y h (x',x 3 ) = Z h (x',hx 3 )-c h . 

(i) There exists y G VF 2,2 (w,M 3 ) such that, up to a subsequence: 

y h —>• y strongly in W 1,2 ^ 1 ,M 3 ). 

The deformation y realizes the compatibility of metrics G and G on the midplate cm 

(3.3) (Vy) T G Vy = G 2x2 . 

The unit normal N to the surface y( oj) and the metric G-induced normal M below have 
the regularity N, M G W 1,2 n L°°(a;,M 3 ): 

\4kt G 


N = 


diy x d 2 y 
\diy x d 2 y\ 


M = 


det G 2 


G~\d iy x d 2 y), 


x2 


where we observe that {diy, GM) = 0 for i = 1,2 and {M , GM) = 1. 
(ii) Up to a subsequence, we have the convergence: 

1 


h 


d 3 y h 


strongly in L 2 (fl 1 ,M 3 ), 


where the Cosserat vector b G W 1,2 n L°°( cc,M 3 ) is given by: 
(3.4) b= (Vy)(G 2x2 )- 1 

(iii) Define the quadratic forms: 


G 13 

G 23 


V det G -> 
+ , _ M. 

det G 2x2 


Q 3 (F) = D 2 W{ld){F,F), 


Q2{x',F 2x2 ) = min {q 3 (g(s)- 1 / 2 #G(s , )- 1/2 ); F G M 3x3 with F 2x2 = F 2x2 } . 
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The form Q 3 is defined for all F G R 3x3 , while Qiix' , •) are defined on F 2 X 2 G M 2x2 . Both 
forms Q 3 and all Q 2 are nonnegative definite and depend only on the symmetric parts of 
their arguments. In particular, when the energy density W is isotropic, i.e.: 

VF e M 3x3 VF e 50(3) W(RF) = W(F), 
then Q. 2 (x', •) is given in terms of the Lame coefficients A, fi > 0 by: 

Q- 2 {x', F 2 X 2 ) = h (O 2 X 2 ) 1 ^ 2 -p 2 x 2 (G , 2 x 2 ) l ^ 2 + ^ + tr ^((? 2 x 2 ) 1 // 2 F 2 x 2 (G 2 x 2 ) 


p2x2 


foT dll -£* 2 x 2 ^ ^sym * 

(iii) We have the lower bound: 


(3.5) liminf ^E h [t?) > T G ${y) := ^ Q 2 (x', (Vy) T G Vb) dx'. 

Proof. The convergences in (i) and (ii) rely on a version of an approximation result from mi; 
there exists matrix fields Q h G W 1 , 2 (o;,M 3x3 ) and a constant C uniform in h, i.e. depending only 
on the geometry of uj and on G, G, such that: 


J |V^(x',x 3 ) - Q'V )! 2 dx<c[h 2 + j^ J dist 2 (G 1/2 V^G- 1/2 ,50(3)) dx), 


J |V</(x ')| 2 dx' < G (l + p J dist 2 (G 1 / 2 Ve /! G- 1 / 2 ,50(3)) dx^j. 

Further ingredients of the proof follow exactly as in [3], so we suppress the details. 


Theorem 3.2. For every compatible immersion y G tF 2 ’ 2 (fl,M 3 ) satisfying (3.3), there exists a 
sequence of recovery deformations G W 1 , 2 (fi /l ,R 3 ), such that: 

(i) The rescaled sequence y h (x',x 3 ) = f h (x',hx 3 ) converges in W 1 , 2 (fi 1 ,R 3 ) to y. 

(ii) One has: 


where the Cosserat vector b in the definition (3.5) of Tq q is derived by (3.f). 


Proof. Let y G TT 2 , 2 (fl,M 3 ) satisfy (3.3). Define b according to (3.4) and let: 


Q = 


d\y d 2 y b 


G w 1,2 n L°°(u, m 3x3 ). 


By Theorem 3.1, it follows that: 

G 1 / 2 QG - 1 / 2 G 50(3) Va.e. x' G w. 
Define the limiting warping field d G L 2 (D,M 3 ): 

" (di b,Gb) 


d{x ) = G~ 1 Q t - 1 c(x', (Vy) T G Vb) - 


(d 2 b,Gb) 

0 


where c(x / ,F 2 X 2 ) denotes the unique minimizer of the problem in: 


VF 2X 2 e K 2 ym Q 2 (x',F 2x2 ) = min|Q 3 (G 1 / 2 (Ff x2 + sym(c <g> e 3 ))G 1/2 ); cGM 3 }. 
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Let {d h } be a approximating sequence in W 1 , 00 (fl,R 3 ), satisfying: 

(3.6) d h -A d strongly in L 2 (w,R 3 ), and h 2 \\d h \\wi,oo ^ 0. 


Note that such sequence can always be derived by reparametrizing (slowing down) a sequence 
of smooth approximations of d. Similiarly, consider the approximations y h £ TT 2 ’°°(w,R 3 ) and 
b h £ W 1,00 {u, R 3 ), with the following properties: 

y h -> y strongly in IT 2 , 2 (w, R 3 ), and b h -A b strongly in VP 1 , 2 (to, R 3 ) 

(3.7) ^ (\\y h \\w^ + ll^lliv 1 - 00 ) < e 

j^\u\ Uh\ —* 0 , where Uh = j x' £ u; y h (x') = y{x') and b h (x') = 6(a/)| 

for some appropriately small e > 0. Existence of approximations with the claimed properties 
follows by partition of unity and truncation arguments, as a special case of the Lusin-type result 
for Sobolev functions (see Proposition 2 in (12]). 

We now define the recovery sequence (f 1 £ W 1,00 (Q h , R 3 ) by: 

i h \x',x 3 ) = y h (x') +x 3 b h {x') + ^-d h (x'). 

Consequently, the rescalings y h £ W 1 , °°(fl 1 ,M 3 ) are: 


y h (x',x 3 ) = y h {x') + hx 3 b h (x') + yxld h (x'), 


and so in view of (3.6) and (3.7), Theorem 3.2 (i) follows directly. The remaining convergence in 
(ii) is achieved via standard calculations exactly as in [3j. We suppress the details. ■ 


It now immediately follows that: 


Corollary 3.3. Existence of a W 2,2 regular immersion satisfying (3.3) is equivalent to the upper 
bound on the energy scaling at minimizers: 

3C> 0 inf E h (£) < Ch 2 . 

£ew 1 - 2 (n h ,R 3 ) 


The following corollary is a standard conclusion of the established T-convergence (see e.g. un¬ 
it indicates that sequences of approximate solutions to the original problem on a thin shell are in 
one-one correspondence to the minimizers of the thin limit variational model X G ( -,. 


Corollary 3.4. Assume that[3.3) admits a W 2,2 -regular solution y. Then any sequence of ap¬ 
proximate minimizers t( h of (3.1), satisfying the property: 


1 


lim yy(ET(n - inf E h ) = 0, 
h ->0 h 2K J 


converges, after the proper rescaling and up to a subsequence (see Theorem 3.1), to a minimizer 
of the functional T G q. In particular, T QG attains its minimum. Conversely, any minimizer y of 

X GG is a limit of approximate minimizers £ h to (3.1). 


On a final note, observe that the defect k(G, G), as defined in (2.8)), is of the order /? 2 min T ( , 
o(h 2 ). It would be hence of interest to discuss the necessary and sufficient conditions for having 
min X G g = 0, which in case of G = M 3 have been precisely derived in [3]. 
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4. AN EQUIVALENT SYSTEM OF PDES FOR (1.3) 


-pm 
mj 1 - iJz-, 


In this section, we investigate the integrability conditions of the system (1.3). Firstly, recall 
that since the Levi-Civita connection is metric-compatible, we have: 

(4.1) dj.G jk = G mk T% 
where the Christoffel symbols (of second kind) of the metric G are: 

r h = + d k G ml - d m G kl ). 

Above, we used the Einstein summation over the repeated upper and lower indices from 1 to n. 

Lemma 4.1. Assume that there is a bilipschitz map f : 12 —> U := f(H), between two open, 
bounded subsets: Cl,U ofW 1 , satisfying (1.3) and such that £ £ W 2 ' 2 (Ct,U). 

Then, denoting by T(j the Christoffel symbols of G, and by T(- the Christoffel symbols of the 
metric G o £ _1 on U, there holds, in Q: 

(4.2) \/i,j, s : 1... n = d m ?fg - o £). 

In particular If automatically enjoys higher regularity: f £ VF 2 ’ 0 O (fl, U). 


Proof. By (4.1) we obtain: 

d<G jk = ((Vfl T G(Vf)) rat f™ + ((V{) T C(VO),„/S, 
while differentiating 0 directly, gives: 

diG jk = ((Vc^fG(VO) + ^ p ((V 0 T 9 p (G o r 1 )^))^ + ((V£) T G(Vc^)) jfc , 

where we used [241 Theorem 2.2.2] to conclude that G o £ _1 £ IE 1 ’ 00 and to apply the chain rule. 
Equating both sides above, and using (4.1) to the Lipschitz metric Go^' 1 , we get: 

d m i s G st d k e f £ + dm^Gstd^T^ 

= dijt'Gstdk? + d^d^dkeiGtmT™ + G sm T™ t ) + d ik f s G s tdj^, 

which we rewrite as: 

0 m £-fg(GV£) sJfc + d m ? f^(GVe) sj 

= dij? {Gvf) sk + djmer ; g {Gvf) sk + d j ?d k z t r* pt (Gvz) 8j + d lk e(GVf) v . 

For each i,j : 1... n we now define the vector Pij £ M n with components: 


(4.3) 


p^ = dij? - d m e f %+dewier. 




■ P9> 


1 


. n. 


By (4.3), we see that: 

Vi, j, k : 1... n 


(Pij, (GVO k—col) — -(Pik, (GVOi- co z), 


where (G\7f) k -col is the vector denoting the fc-th column of the matrix GV£. 


Since P^ = Pji, it now follows that: 


(4.4) Vi, j,k : 1 .. .n (Pij, (GVOk-col) = 0 . 

But the columns of the invertible GV£ form a linearly independent system, so there must be 
Pij = 0, which completes the proof of (4.2). ■ 
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Corollary 4.2. Let £ : fl —>■ U be a bilipschitz map satisfying (1.3) as in Lemma f.l Assume 
that £ E VL 2,2 (17, U). Then £ and £ = £ _1 are both smooth and bounded, together with all their 


derivatives. In particular, (4-2) holds everywhere in 17. 

Proof. It was already established that V 2 £ E L°°. We have F := V£ = (V£ o £) _1 E L°°. Notice 
that F _1 is the composition of two Lipschitz mappings and hence it is Lipschitz. We conclude that 
for all i, djF = —F8i(F~ 1 )F E L°°, which implies that VC E W 1,0 °, and hence C € W 2 ,ao (U, 17). 
By (1.3) we get the following formula: (V£) T (G o C)(VC) = Go (, valid in U. By the same 


calculations as in Lemma 4.1, it results in: 


Vi, j, s : 1 ... n diiC = d m m - ckf'OjCn- 


s 

vr 


In view of (4.2) and by a bootstrap argument, we obtain that £, C E W k,ao for every k > 1. Hence 


the result follows. 


Corollary 4.3. Equivalently, (4.2) can be written as: 


Vi,j,s:l...n 8^ = 8^% 

+ 8 j ed l G mp - dj^di^dtGpqilyi)- 1 ) 


tm 


Proof. Denoting C = £ 1 as before, we obtain: 


ps 1 r~ism 

1 m ~ 2 ^ 


dpCdtGmq + dqCf 8tG m p — dmCf'df Gpq^ . 


Inserting in (4.2) and contracting dpCfdjf? to the Kronecker delta ST we obtain (4.5). 


Theorem 4.4. Consider the following system of the algebraic-differential equations in the un¬ 
knowns £, : 17 —>■ M n , i : 1... n: 

Gij = wjGstWj 
w\ = diff 


(4.6a) 

(4.6b) 

(4.6c) 


djW? = '<f- - \G sm [w\d 3 G rnv + wfjdiGmp - W^dtGpqW, 


t 

m I i 


where the matrix \W( n \t,m=i..n is defined as the inverse of the matrix field w : ft —>• M nxn , whose 


columns are the vectors Wi, i.e.: Wf = ( Ya =i w% ® gY™■ Then we have the following: 

(i) Problem 0 has a solution given by a bilipschitz map C E W 2 , 2 (Q,U) (as in Corollary 
1 4-2j ), if and only if (4.6) has a solution (£,w) given by a bilipschitz £ : 17 -A U and 
w E W 1,2 (17,M nxn ). 


(ii) Problem ( |4.6| ) has a solution, understood as in (i) above, if and only if (4.6c) has a 
solution. This statement should be understood in the following sense. 


Assume that (4- 6 c) is solved in the sense of distributions, by the vector fields Wi E 
L°°(17,M n ), i : 1.. .n, such that W ^ are well defined and Wf E L°°(17). Then there exists 
a smooth £ : 17 — > U such that (4- 6 b) holds. Moreover, £ is locally invertible to a smooth 


vector field (, and the Christoffel symbols of the following metrics: 

G = Go( and g = (V() t (Go()(V() 
are the same. If additionally £ is globally invertible to £ : U —>• 17, and if we have: 
(4.7) 3 xq E 17 G(xo) = uj t Gw(x o), 















THE METRIC-RESTRICTED INVERSE DESIGN PROBLEM 


13 


then (4-6a) holds in 12 . 


Proof. Clearly, the equivalence in (i) follows from Lemma 4.1 and Lemma 4.3 For the equivalence 
in (ii) note first that, by a bootstrap argument, an L°° vector field w satisfying (4.6c) is automat¬ 
ically smooth in 17, together with W. Further, in a simply connected domain 17, the condition 
(4.6b) is the same as: 


Vi, j,t : 1.. .n 


djw\ = diwj, 


which is implied by (4.6c), by the symmetry of its right hand side (in i,j). 

We now prove that the metrics G and G have the same Christoffel symbols on the subdomain 
of U where the local inverse £ is defined. Note that both £ and £ are smooth. We first compute: 


(4.8) 


dm di( P G ps dj( S = dimCGpsdjC + diCGpsdjmC + djCd m Cd q G 


'qvps- 


By Lemma 4.1, we also obtain: 


djwf = w s m fg - w p j w q i T s pq , 


where pL are the Christoffel symbols of Q. Now, since: d p Q s w p = Sf = diCfWp, it follows that: 

d pq Cwfw p = -dpCdjw f = -d p C ('<f" l - wfw?T p af) } = -f £ + d p ( s w?w?r p afj . 
Consequently, (4.8) becomes: 

dmGij = dimCGpsdjC + diCGpsdjmC 


+ d^djCdm^Gapf ~ d^CwZw* + dtCw^T^ 


(4.9) 


+ 


diC p djCdmC q Gas(^- d^sCvJqWp + dtCw^WpT^ 


= diCdtCGapTInj + djCdtCGasT 
= GtiP m j + GuT. 


t 

im 


3 mi’ 


Call 7 - the Christoffel symbols of G, so that: d m Gij = GtiJmj + Gtjlmi- Therefore: 

Vi, j, m : 1 ... n Gtjlmi ~ 7L) = -GtiiJ'mi ~ itni)- 


Using the same reasoning as in the proof of (4.4) we get, as claimed: 

~\m 
■ ji 


Vi, j,m : 1 ... n 


yrn _ m 
X 11 111 ’ 


This concludes the proof of (ii) in view of Lemma 4.5 below, and since G(£(x o)) = £/(£(zo))- 


Lemma 4.5. Assume that the Christoffel symbols of two smooth metrics Gi,Gi on a connected 
domain U C M n coincide, and that for some xq E U, Gi{xq) = G 2 (xo). Then G 1 = G 2 - 

Proof. Let V* represent the covariant derivative associated with the metric G% through the Levi- 
Civita connection m p. 114 and Theorem 2.2, p. 158]. An immediate consequence of the 
equivalence of the Christoffel symbols (in the identity coordinate system of U) is that V := Vi = 
V 2 [13 p. 144], By |T3 Proposition 2.1, p. 158], we have ViGi = 0, i.e. both Gi and G 2 are 
parallel tensor fields on U with respect to V m pages 88 and 124], Let 16 !/ and let 7 be any 
piece-wise C 1 curve connecting xq and x in U. Let r 7 be the parallel transport from the point uq 
to u along 7 . It follows that: 


r 7 (£i(x 0 )) = Gi{x) Vi = 1, 2, 
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which concludes the proof. ■ 

Remark 4.6. Assume that the Christoffel symbols of two smooth 2d metrics Q\ . Q 2 on a connected 
domain U C M 2 coincide. Assume moreover that the Gaussian curvature of Q\ (and hence of Q2) 
does not vanish identically in U. Then, one can directly prove that there exists a positive constant 
A such that Q\ = \Q-2 in U. 

Namely, let Vj represent the covariant derivative associated with the metric Qi through the 
Levi-Civita connection, as in the proof of Lemma |4.5| Again, equivalence of the Christoffel 
symbols implies that V := Vi = V 2 . Let xo G U and let H = Hol(17, xq) be the holonomy group 
associated with V at xq. By [5], (see also [21 Theorem 392]), H is a connected subgroup of the 
special orthogonal group associated with the scalar product Qi{x 0 ). Since (H, Qi) is not flat, H 
cannot be the trivial subgroup. The only other possibility (when n = 2) is that H is the entire 
50(2, Gi(xo)). This implies that the SO( 2 ,Qi(xo)) = 50(2, Q 2 (xq)), which by the transitivity of 
the action of the special orthogonal group over the unit sphere, results in Qi(xq) = Af^Oco)- 

Note that for n > 2, the holonomy group H generically coincides with the full special orthogonal 
group [21 p.643]. In this case again SO(n,Q i(xo)) = 50(n, £ 2 (^ 0 )) for some xq G U, and so the 
result is established for generic metrics Qi. 


5. Some remarks and examples 


Remark 5.1. Note that in problem (1.3) one can, without loss of generality, assume that: 


(5.1) G(x 0 ) = G(x 0 ) = Id n 

Indeed, denote A = G(x o)~ 1//2 and A = G(x o)~ 1//2 the inverses of the unique symmetric square 
roots of the metrics at xo, and let Hi = A _1 (H) be the preimage of H under the linear transfor¬ 
mation x 1 —>■ Ax. Then G± = A(G o A)A and G 2 = A(G o A)A are two metrics on Hi, that equal 

W 


Id at x\ = A 1 xq. Further, (]l.3h can be written as: 


(5.2) Gi(x) = (a-^v^Ax))!)^!^)^- 1 ^^))!) = V6GiV6(x), 


where ^1 = A (o A. Clearly, existence of a solution ^1 to (5.2) on Hi is completely equivalent 


with existence of £ solving (1.3) on H, with the same required regularity. 


Notice also that, in view of the compensated regularity for £ in Corollary |4.2[ any solution to 


(1.3) satisfies, up to a global reflection, the orientation preserving condition (H.lk In view of (5.1), 


(4.7) now becomes: w(x 0 ) G SO{n). For the particular case of w(x 0 ) = Id, it is easy to notice 
that (4.6c) at xo reduces to: djwf(x 0 ) = Tfj — xf 
the metric G on H. 


Ur 


with xij denoting the Christoffel symbols of 


Remark 5.2. We now present alternative calculations leading to the system (4.6a) - (4.6c), adapt¬ 
ing classical ideas 
of the connection. 


These calculations do not require knowledge of the metric compatibility 


1. Assume that the re ex ists a bilipschitz map £ : H —> U, whose global inverse we denote by 
£ : U —> H, such that (1.3) holds in H. Let F^fc the Christoffel symbols of the first kind of the 
metric G onH: 

1 'ijk o (J-^iGjk T djGik d k Gij), 
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while let Y a p^ stand for the Christoffel symbols of the metric Q = G o ( on U: 

Y aP 7 — 2 (^aGp'y "F dpGa-y d'yGap') • 

As above, w e wi ll use the Latin indices for components of vectors in f2, and Greek indices in U. 
Since by (1.3) we have G = (V£) T (f? o £)(V£), we obtain: 

d k G ij = d k {d i i a d j ^{G a p oO) 

= (dkiCdj^ + d kj Cdi^ (Gap o £) + diCdj^ (( d^Gap ) o ^ d k C- 


Similarly, exchanging the following pairs of indices: 7 with a and (3 with 7 in the first equation, 
and 7 with /3 in the second equation below, we get: 


diG jk = (dijCdk^ + d ik CdjC 13 ^ (Gap ° 0 + d^ a dj^ {(d a G 1 p) o ^jd k C, 
djG ik = (dijCdk^ + d jk Cdi^ (Gap ° 0 + diCdj^ {(dpG ai ) ° £ )d k C. 
Consequently: 

(5.3) f = d^djtfdk? (r q/37 o e) + 4 ^ (^Gap o e). 


2. We now compute the Christoffel symbols of the second kind for the metric G: 

■p/c pfem-p 

-L ^ ^ ijm- 


Since G 1 = ((VC)<7 1 (VC) T ) 0 £, it follows that: G km = (d a ( k dp( m G a ^) 0 £• By (5.3) we get: 

(5.4) 

Y% = [(d P c k d u c m G pu ) o c) o o + ((a P c*a„c m £n o e) (G a p o 

= ((£ PT r a p 7 ) o tyiFdjt^dp c k o 0 + ((^ aj8 ) o (a P c fc o 0 

= ( r a/3 0 medj^(d p C k o 0 + (daC k o Odijt a , 

where we contracted <9 m £ 7 (d v C, m o £) to the Kronecker delta SI and d m ^(d u ( rn o £) to 8 ® in the 
second equality, and G p/3 G a p to 5a in the third equality, where we also used the definition of the 
Christoffel symbols r p ^ of the second kind of the metric G■ Further: 

t%d k e = (OaC k o medijc + (r^ o oa^a^^c* ° o 

= ^ + (^0^3/, 

which implies the same formula as in (4.2): 

(5.5) d v e = V'l/h-e - Kp o Odi^dj^, 

3. We now proceed as in Corollary |4.3| 

r «/3 = (daGp-, + dpGa~t - Gap} 

= \G^[(d s Gp 1 O CRC S + (drGaj o C) dpC - (d q G a p o c)a 7 c 9 ) . 
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Consequently, and in view of (5.5): 


f %d k e - = (Kp o 


= -{gn o 0 [d s G^{d a Q s o mrdjf + d r G ai (d 0 C o OdiCdj^ 

- d q c a ^c q o medjt 0 

= \g^ (dsGftSidjtf + drG^diC - dgGapW^diCdj^) , 


which directly implies (4.6c). 


Remar k 5 .3. We further observe that if (4.6c) hold s, then defining £ by (4.6b) we may again 
obtain (5.4) simply by reversing steps in R emar k 5.2 Letting G = (V£) T (<7 o £)(V£) and going 


through the same calculations, the formula (5.4) follows, this time for the Christoffel symbols L'L 
of G, i.e.: 


f& = (r^ o mcdj^(d p c k o 0 + ( d a ( k o Qdij? 


4.5 


Hence, we see that T - = Tfj in H and again, in view of Lemma 

(V£) t G(V£) in H, provided that we have this identity at a given point xq, as required in (4.7). 


it follows that G = G = 


Example 5.4. Assume that G is constant. Then the equations in (4.6c) become: 

(5.6) djwf = w s m f%, 

whereas the Thomas condition (eh, see also next section) for the above system of differential 
equations is: 

<4f Tj + vfilrrFTj = w m d ^Tk + w a p r p jm r ? k Vi,j,s,k: 1 

Equivalently, the following should hold: 

fcfE - dif s + f zf - rzrl) = o v», j, s, k ■. 1 ... n , 


, n. 


UK 


(^k 1 ij (; j 1 ik L kp L ij L jp^ ij 


at all points x in a neighborhood of x$, and for all w in a neighborhood of a given invertible 
wq £ M nxn which satisfies G(x o) = WqGwq. Since the expression in parentheses above equals 
K!jjk( x ), it follows that the Thomas condition for (5.6) is precisely the vanishing of the whole 
Riemann curvature tensor R... of the metric G, or equivalently that G be immersible in M n . 

On the other hand, letting A = \[G denote the unique positive definite symmetric square root 


of the matrix G, we see that the problem (1.3) becomes: 

G = (V£) r GV£ = (AV£) r (AVO = (V6 ) T V£i, 


with £i(x) = A£(x). Since A is invertible, we easily deduce that (1.3) has a solution iff G is 
immersible in M n . Consequently, in this example the Thomas condition for (4.6c) is necessary and 
sufficient for solvability of (El- 


Example 5.5. Assume that G{x) = /t(x)Id 3 . Then the problem (1.3) becomes: 


(5.7) 


(VO T V£ = 1 G = e 2f G, 


with / = - -log//. 


Solution to (5.7) exists if and only if the Ricci curvature of the metric e 2 ^G, which is con¬ 


formally equivalent to the metric G, is equal to 0 in fL More precisely, denoting: V%f = 

Cr 
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3x3 , A Q f = & k dj k f - & k T l jk d l f G R, and ||V/||| = & k djfd k f G 


A,f - r %d k f] iJ:1 .. 3 G 
the condition reads: 

(5.8) 0 = Ric(e 2 ^G) = Ric(G) - (V|/ - V/ ® V/) - (a^/ + || V/|||) G. 


When also G = A(x)Id 3 , then (5.8) after setting h = \ log reduces to: 


(5.9) 0 = -(V 2 /i- Vh® Vh) - (A/i + ||V/i|| 2 )ld 3 . 

An immediate calculation shows that the only solutions of ( |5.9[ ) are: 

h(x) = —21og|x — a\ + c or h = c, 


with arbitrary constants c G M and a G M n . Indeed, in this case the solutions to (5.7) are 
conformal, and so by Liouville’s theorem they are given by the Mobius transformations in M 3 , as 
compositions of rotations, dilations, inversions and translations of the form: 

R(x — a) 


£(x) = b + a- 


\x — a\P ’ 


R G 50(3), a, 6 G M , a G M, /3 G {0, 2}. 


Example 5.6. In dimension n = 2, existence of a solution to (5.7) is equivalent to the vanishing 
of the Gauss curvature: 


(5.10) 


0 = K(e 2f G) = e 2 / (k(0) - A^/), 


which further becomes: k(G) = A gf. When also G = A(x)Id 3 , this reduces to: Ah = 0, which is 
also equivalent to the following compatibility of the Gauss curvatures of G and G: = j. 

6. Some further remarks on systems of total differential equations 

Systems of total differential equations have been extensively studied in the literature mm® 
mm- They are over-determined systems in the unknown w : —> M . N , of the form: 


( 6 . 1 ) 


diw(x) = fi(x,w(x)), 


Vi = 1 ... n Vx G O, 


where all the partial derivatives of w are given by functions fi : 12 x 


Note that if 


the latter are assumed of sufficient regularity, the uniqueness of solutions to (6.1) with a given 


initial data w(x o) = wq is immediate. For existence, observe that any solution must satisfy the 
compatibility conditions: 

dijW = djiiv. 

This leads, under sufficient regularity assumptions, to the necessary condition for existence of w: 

( 6 . 2 ) 


— + = Vi, j : 1... n VxGh 


V dxj dxi ' v dw' 

The advantage of the system ( |6.2[ ) is that it does not involve any partial derivatives of the a-priori 
unknown solution w, and hence, if certain conditions are satisfied, it can be used to obtain the 
candidates for w by solving for w(x) at each x G fi. 


Naturally, the more solutions (6.2) has, the more there is a chance to find a solution to (6.1). A 


plausible strategy, possibly adaptable as practical numerical schemes, is to find all the candidates 


w from (6.2) and check whether they satisfy (6.1). This insight, combined with the observation 
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about the over-determination of the original system, implies its rigidity, and leads to non-existence 
of solutions in generic situations. On the other hand, the ideal situation is to have: 


(6.3) 


dfi dfj ,dfi s /dfj \ , 

w s - + («£)/> 


dxi v dw J v dw' 
satisfied for all x E 0 and all w E M^. All functions re : hi 


Vi, j : 1. .. n. 

obtained this way are solutions 
to ( 6 . 2 ) and, as shown in | 21 j . this leads to existence of a family of local solutions to ( 6 . 1 ), 
parametrized by the initial values at a given point xq E 0. If (6.3) is not satisfied, one comes 
short of having an ample set of solutions w and might be content for other intermediate scenarios, 
where solutions to ( 6 . 1 ) may still exist but with less liberty in choosing initial values. 

In section [4] we showed that (1.3) is equivalent to a system of total differential equations. 
Note that even though our problem shares some familiar features with the isometric immersion 
problem, it is of a fundamentally different nature in as much as we cannot establish the equivalence 
of (6.1) and (6.3). Not being able to close the system (6.2) as in the isometric immersion case, in 
order to find necessary and sufficient conditions, we need to study the above mentioned possible 
intermediate scenarios when the Thomas condition (6.3) is not satisfied. We will carry out this 
plan in this sections [7] and [8j Below, we begin by a simple example of (6.1), whose conditions 

(6.3) are far from optimal. 

Example 6.1. For n = 2,1V = 1, consider the system: 

(6.4) Vtc = w 2 a + wb + c in 12 C M 2 , 

where a, b : 12 -A M 2 are given smooth vector fields and c E M 2 . In order to find the Thomas 
condition (6.3) for existence of a solution w : 12 —> M, we take: 

0 = curl(u> 2 a + wb + c) = rc 2 curl a + 2w(V' L w, a) + rccurl b + (V^ru, b). 


Substituting (6.4), we obtain the counterpart of (6.2) in the present case: 

(6.5) (curl a + (a, b^^w 2 + (curl b)w = 0. 


The satisfaction of the above in the (x,w) space is precisely the condition (6.3): 

( 6 . 6 ) curl a + (a, h 1 ) = 0 and curl b = 0 in 12 . 


Lemma 6.2. Condition (6.6) is not necessary for the existence of solutions for ( 6 . 4 ). Also, 
existence of a solution to (6.5) does not guarantee the existence of a solution to (6.f). 


Proof. 1. Let w : 12 —> M be any positive smooth function. Let c = 0 and let b : 12 
vector field which is not curl free. Define: 

a = —^(Vmi — wb). 
vr 


be a 


Then ( |6.4[ ) is satisfied but not ( | 6 . 6 [ ). 

2. Let a, b : 12 —> M 2 be two smooth vector fields such that the first condition in ( 6 . 6 ) does not 
hold at any point x E 12. For example, one may take: 

a(x) = x and b(x) = x^~, 


in any domain 12 which avoids 0 in its closure. Then (6.5) becomes: 
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and it has only two smooth solutions: w = 0 and w(x) = —2/|x| 2 . On the other hand: 


V(-, 


\x\ 


. 2 ,2 2 j_ 

F "(ul) X ~ 02 X + O 


so none of these functions is a solution to (6.4) when c / 0. 


Further, observe that augmenting (6.4) to a system of (decoupled) total differential equations: 


Vic/ = w'jdj 4- wjbj + cj, 


I = 1... N, 


(6.7) 

gives the condition: 

(6.8) (curl ai + (a/, bj~))w‘j + (curl bj)iuj = 0 V/ = 1... N, 
resulting in: 

(6.9) curl aj + (a/, bj~) = 0 and curl bj = 0, in 14. 


As in Lemma 6.2 we can set up the data such that, there exists a solution w±(x) to the first 


equation in the system (6.8) and that all the equations in (6.9) for I = 2 ... N are satisfied in such 


a manner that (6.7) still has no solutions. 


7. The Thomas condition in the 2-dimensional case of (1.3) 


In this section, we assume that n = 2 and follow the approach of [7j. Recalling (2.3) as the 
equivalent form of (1.3), we note that it has a solution (on a simply connected 14) if and only if: 

(7.1) curl(G- 1/2 RG 1/2 ) = 0 in 14, 

for some rotation valued field R : 14 -» 50(2). Denote V = G 1 / 2 , V = G 1 / 2 and represent R by a 
function 8 : 14 —>• M, so that R(x) = R(8(x)) is the rotation of angle 8 . Note that: 


7r, 


djR{8) = (dj8)RW, where W = R{~) = 


0 -1 

1 0 


Since IT is a rotation, it commutes with all R E 50(2). Also: W 1 = W 1 = —W, and: 

(7.2) VF E GL(2) WFW = —cof F = — (det F)F~ l . 

We finally need to recall the conformal-anticonformal decomposition of 2 x 2 matrices. Let 
M 2x2 and M 2x2 denote, respectively, the spaces of conformal and anticonformal matrices: 


»2x2 


It is easy to see that 


a b 
—b a 
2x2 _ in 


; a, b E 


»2x2 


a b 
b —a 


a,b E 


2x2 


d 2x2 because both spaces have dimension 2 and they are 


mutually orthogonal: A : B = 0 for all A E 


For F = [Fij]ij : i ,2 E M 2 , its projections F c on 


v ^ x2 and B E M 2x2 
2x2 


, and F a on are: 


2x2 


1 


F c = - 
2 


Fn + F22 F \2 — F 2 i 
F 21 — F\ 2 Fn + F 2 2 

We easily obtain that: 

(7.3) F = F C + F a , |F| 2 = | F c | 2 + |F a | 2 
Also: 

(7.4) VF E 50(2) 


F a = - 
2 


Fn — F 2 2 F\2 + F 2 i 
F12 + F 2 i F 2 2 — Fn 


and detF = 2(|F C | 2 — |F a | 2 ). 


F C R = RF C and R T F a = F a R. 
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Lemma 7.1. The following system of total differential equations in 6 is equivalent with (1.3): 


(7.5) V0 = —— 
det V 

where we have denoted: 


- V curl V + 


1 


det V 


V A 1 Ve 2 -A 2 Ve 1 ) + 


1 


(7.6) 


Ai = ( VdiV - 1 ) c 


det V 

B i = (VdiV~ 1 ) a Vi = 1,2. 


VR(-29) [B l Ve 2 -B 2 Ve x , 


Proof. We calculate the expression in (7.1): 

curl (V-'RV) = d^V-'RVef) - d 2 {V~ 1 RVe 1 ) 


= (9i V-'RV + V~ 1 Rd 1 V)e 2 - (, d 2 V~ l RV + 
+ (SiejV'-^WV'ea - {d 2 9)V~ l RWVei 


The two last terms above equal, in view of (7.2): 

(d l e)V~ 1 RWVWe l + {d 2 9)V~ 1 RWVWe 2 = ( V~ 1 RWVW)(y9) T = (det V) (V- 1 RV~ 1 ){V9) T . 
Therefore, by ( |7.1[ ) and ( |7.4[ ): 

V9 = -j^(VR t V) ( - (5i V-'RV + V~ l Rd 1 V)e 2 + {d 2 V~ 1 RV + V~ l Rd 2 V)e i) 


1 


VR 1 Vd 1 V- 1 RVe 2 - V- l R 1 Vd 2 V- L RVe 1 + V~ L diVe 2 - V~ l d 2 Ve i 


r— 1 dT 


-1 : 


r—1 1 


r—Xi 


det y 

We simplify the equation for 9 as follows: 
1 


V0 = 


det y 


VA 1 ye 2 - yA 2 y ei + VR(-29)B 1 Ve 2 - VR(-29)B 2 V ei 
+ Vd 1 Ve 2 -Vd 2 Ve 1 ), 


which clearly implies (7.5). ■ 

If the Thomas condition for ( |7.5| ) is satisfied, there exists a unique solution to (7.5) for all initial 
values 9(x o) = 9q. The original problem B has then a unique solution for all initial values of 
the form w(xq) = V£(xq) = G~ 1 / 2 RqG 1 ^ 2 (xo) with Rq £ 50(2). 


Example 7.2. Assume that G = M 2 . Then A{ = Bi = 0, so that (7.5) becomes: 

1 


(7.7) 


V9 = h 


v 


where hy = 


det V 


-VcuvlV. 


The condition for solvability of ( |7.7| ) is: curl hy = 0, which is expected from Example 5.4 Indeed, 
a direct calculation shows that (see also [SJ Remark 2, page 113]) this condition is equivalent to 
the vanishing of the Gaussian curvature of G: 

(7.8) k(G) = -—curl hi. 


det V 




Example 7.3. As in Examples 
(O). Then: 


5.5 


(7.9) 


and 


5.6 


assume that G = qi(x)Id 2 = e 2 ^Id 2 , with / given 


m 


1 


VdiV = = (dif) Id 2 . 
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This implies that Ai = (<9j/)Id 2 and Bi = 0. Consequently, (7.5) becomes: 

VO = —Lycuriy + —^-^t^ 2 ((o> 1 /)e 2 - (d 2 f) ei ) = hy + ~jL=G V ± f. 
det V det V V det G 


But by (7.2) it follows that: 

1 


\J det G 


G V ± f = . 1 _ GWVf = 1 _ W T WGWVf = V det G WG^Vf 


\J det G V det G 

(VdetG G -1 Vf) ± , 


and hence: 


curl ( 1 G V ± f) = div('\/det G G _1 V/) = VdetGA^f. 

V ydet G 


Thus, the Thomas condition to (1.3) is here: 


curl hy + V det G A yf = 0, 


and we see that, in view of (7.8), it coincides with the equivalent condition (5.10) for existence of 
solutions to (1.3). 


Lemma 7.4. Without loss of generality and through a change of variable, we can assume that: 


G = A(s)Id 2 = e 29 Id 2 


with 


9 = ^logA. 


Then, for an arbitrary metric G, condition (7.5) becomes: 
(7.10) VO = m + R(-20)n, 


where m = V^g + (A\e 2 — A 2 e\) and n = Bie 2 — B 2 ei. The Thomas condition of (7.10) reads: 
(7.11) 


curl m — 21 n | 2 = 0, 
divtt - 2(rr L ,m) = 0 , 
curl it — 2(tt, m) = 0. 


Proof. As in (7.9), we observe that VcurlV = e 29 V~ L g, and hence: 

V0 = V^g + (Aie 2 — A 2 ei) + R(—2Q) ( Bie 2 — B 2 e\). 

The above equation has a similar structure to (24) in [7], even though the two original problems 
are different. In the present case, both G and G are involved in defining m, while in [7] the vector 
fields m, n depend on the matrix field G in the Left Cauchy-Green equation: ( Vr])(Vr]) T = G. 

In order to derive the Thomas condition for (7.10), note that: 


curl V 0 = curlm + di(R(—20)n, e 2 ) — d 2 (R(— 20)n, e\) 

= curlm + (R(—20)din, e 2 ) - {R(-20)d 2 n, ei) - 2d 1 0(R(-20)Wn, e 2 ) + 2d 2 0{R{-20)Wn, ef). 





























22 


AMIT ACHARYA, MARTA LEWICKA AND MOHAMMAD REZA PAKZAD 


Substituting 8$ from ( |7.10 ) we arrive at: 

0 = curlm + {R(—26)d\n, e 2 ) — {R(—29)d2n, ei) 

- 2(m + R(-20)n, ei){R{-28)Wn, Wei) - 2(m + R(-29)n, e 2 )(R{-29)Wn, We 2 ) 
= curlm + (R(—20)din, e 2 ) — (R(—28)d 2 n, e\) — 2{m + R(—20)n, R(—29)n) 

= curlm + (5m, R{26)Wei) + (5 2 n, R{26)We 2 ) - 2|n| 2 - 2(m, R(-28)n) 

= curlm - 2|n| 2 + (Vn : R(29)W) - 2(m, R{-20)n) 

= curlm — 2|n| 2 + (Vn — 2Wn ® m : R(29)W). 

Finally, writing R(29)W = R(28 + ir/2) = — sin(20)Id2 + cos(2 9)W, we obtain: 

0 = curlm — 2|n| 2 — sin(20)(divn — 2(tr L ,m)) + cos(20)(curln — 2(n,m}), 


which should be satisfied for all x G and all 9 G [0, 2ir), implying hence (7.10) 


Example 7.5. In the setting and using notation of Example 5.6, we see that: m = V^g + V x / = 
V _L (/ + g) = V ± h where h = ^ log j). and tt = 0. The Thomas condition for (1.3) here is hence: 
curlm = 0, that is equivalent to: 

Ah = Alog( —) = 0, 

9 


which is further exactly equivalent to existence of solutions in Example 5.6 


Example 7.6. We will n ow p rovide an example, where the Thomas condition (7.11) is not sat¬ 
isfied, but a solution to (1.3) exists. We start with requesting that G = M 2 , which yields: 


m = A\e 2 — ^ 2 ei and n = B\e 2 — B 2 e\. Consider the general form of the diagonal metric G : 

e 2 a{x) q 

0 e 2b ^ 

for some smooth functions a, b : £1 —> M. Then: VdiV~ l = — diag{<9ja> 5*6}, and hence: 


G{x) = 


which leads to: 


At = ~^di(a + fr)Id 2 , 


m = —-V _L (a + b ) 


Bi = ^diag {di(b - a), di(a - b )}, 


and n = ~(d 2 (a - b), 5i(a — b)). 


Therefore, the Thomas condition (7.11) becomes: 


A(a + 6) + |V(a — b)\ 2 = 0, 

5i2(a - b) + \di(a - b)d 2 (a + b) + \d 2 (a - b)d±(a + b) = 0, 

5n(a - b) - d 22 (a - b) + 5i(a - b)d\(a + b) - d 2 (a - b)d 2 (a + b) = 0. 

It can be simplified to the following form, symmetric in a and b\ 

' A(a + 6) + |V(a — b)\ 2 = 0, 

(7.12) < 25i20 + 5ia52a = 2d\ 2 b + dibd 2 b, 

5n a - d 22 a + (5ia) 2 - (5 2 a) 2 = dub - d 22 b + (5i6) 2 - ( d 2 b ) 2 . 
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We now specify to the claimed example. Take a = — logxi and b = — logx 2 , so that: 

-2 


G( Xl ,x,)= J x 


0 

-2 

2 


defined on a domain ff C M 2 in the positive quadrant, whose closure avoids 0. We easily check 
that the first condition in (7.12) does not hold, since: 

Alog(xix 2 ) + |Vlog( —)| 2 = -d n logxi + |di logxi| 2 - <9 22 logx 2 + |<9 2 logx 2 | 2 = 2(\ + \). 


X2 




On the other hand, Co(xi,X 2 ) = ^(xf,x 2 ) solves (1.3), because: 
(V(o) t GV( 0 = 


1 

o 

H 

i_ 


x 1 2 0 


l- 

O 

i—i 

0 x 2 


0 X 2 2 


0 x 2 


= Id 2 = G. 


In fact, by reversing the calculations in the proof of Lemma 7.4 in can be checked that £o and 
—£o are the only two solutions to the problem (1.3). 


8. A SUFFICIENT CONDITION FOR THE SOLVABILITY OF THE PROBLEM (4.6C) (4.7) 


In this section we go back to the general setting of n > 2 dimensional problem (1.3) and its 
equivalent formulation (4.6c). Let wq G M nxn be an invertible matrix, such that: 

G(x 0 ) = ivqG(xq)w 0 . 

For every s,i,j : 1... n let f s f‘ : O —> M be the smooth functions of (x,w), dehned in a small 
neighborhood O of (xq,wq), and coinciding with the right hand side of (4.6c): 

fj’ 2 (x, w) = w s m f™ - ^G sm (w r -d J G mp + w p -diG mp - 


Jm I » 


where in” 1 denotes the element in the s-th row and m-th column of the matrix w. 

Differentiate fj' 1 formally in d/dxk, treating w as a function of x, and substitute each partial 
derivative d}~Wq by the expression in /£’ 9 . We call this new function and note that for w 
satisfying (4.6c) one has: dkdjwf(x) = Fj%(x,w(x)) for every x G D. 

Now, Thomas’ sufficient condition [21] for the local solvability of the total differential equation 


(4.6c) with the initial condition (4.7) where w{x o) = wq, is that: 

, n 


( 8 . 1 ) 


Vs, i, k,j : 1 


Fj£ = F°£ - F%’j = 0 in O. 


Following [I], we now derive another sufficient local solvability condition. For convenience of the 
reader, we present the whole argument as in section 4 [1]. 


For every a\ = 1 .. .n, let = 0 be the equation obtained after formally differentiating 


F 


= 0 in d/dx ai and replacing each partial derivative d ai wj by fa\ as before. Inductively, for 
every m > 0 and all a \, a 2 ... a m : 1... n, we define the equations Ff’l „ „ =0. 

...LX-m 

Theorem 8.1. Fix N > 0. Assume that there exists a set S C {1... n} 2 of cardinality 1 < #S = 
M < n 2 , and there exist M equations F 1 = 0, F 2 = 0, ... F AI = 0, among the equations: 

(8.2) {j# = o}. iWit u { f% m = u... u =o} s . iAai 


a N ' 
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such that: 
(8.3) 


det 


dF r 

dWq 


(xo,w Q ) 


(p,q)£S 


/o. 


From now on, we will denote by w p the coefficients in the given matri x w with (p, q) E S and by 
w[ the coefficients with ( l,t ) 0 S. By the implicit function theorem, (8.3) guarantees existence 
of smooth functions v p : U x V — > M for every ( p,q) E S, where {v p (x, {w L t\ (p,q)es G 
is defined for x in a small neighborhood U of xq and for { v [} in a small neighborhood V of 
{Oo)l}(z,t)£S> satisfying: 

V(p, ?)G5 vjj(x 0 , {(tD 0 )i}) = (w 0 ) p , 

(8.4) 

Vr : 1... M VxeU V{vj} E V 
Assume further that: 

\/m:0...N + l \/i,s,j,ai...a m Mx&U V{uJ} 6 V 


{w?(*,{vt}(i J t) 0 s)} (Pi , )eS ) = 0. 


(8.5) 


F 


j,k,a 


x i{ v t},{ V q( X ’i v t 


= 0. 


Then the problem (f.6c) (f.l) has a solution, defined in some small neighborhood U of xq, such 
that w{x o) = Wq. 

Proof. 1. Below, we drop the Einstein summation convention and use the sign instead, for 
more clarity. By U and V we always denote appropriately small neighborhoods of xq E M n and 
{(^o)t}(z,t)0S G , respectively, although the sets U and V may vary from step to step. 

We will seek for a solution in the form ^{w)|(x)}, {wq(x) = v p (x, {u)|(x)})}^ , where the functions 
v p are defined in the statement of the theorem. Note first that by (8.5), for every a \ 1... n, every 
m : 0 ... N, and s, i, j, k, a\ ... a m , after denoting: F = Fj'(. Qi we have: 

r)TT „ AT? . . 

(8.6) 


^771 _ _ 1377' 

(x,u)+ ^l( x i v )fa( x i v ) + Y A T3>( X ’ V )fZ q ( X i V ) = 0 ’ 


dx 


( 1 M s 9t6 ‘ 


(p,q)es 


for all x E U and all {v l t } ( i^ s G V, where above v = {{v l t }(ij)£s,{ v q}( P ,q)es} and Vq(x) = 
Vq(x,{v l t (x)}) for all (p,q) E S. 

2. Fix now any point x E U and let [0, 77 ] 3 r i-a g(r) = {g a (r)} a= \,, /n <3 U be a smooth 
path with constant speed, such that: < 7 ( 0 ) = x. Consider the solution [0,??] 3 t i-a h(r) = 
G V of the following initial value problem: 


(8.7) 


dT 

3 


-TZ h t( T )= Y fa(9(T),h(T),{vP(g(T),h{T))}Y^-(T) Vr E [0,1] V(/,f)^5 


ot=l...n 


dr 


. K(0) = 4- 


Applying ( |8.5[ ) , we easily see that for F denoting, as in step 1, any function in the set of equation 
(| 8 . 2 l), there holds: 


F 


(ff(r),h(T), {v%(g(r),h(T))}) =0 Vr E [0,rj]. 
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Differentiating in r yields: 

dl 

dx a 


a=l...n 


dr 1 (I®** 


E 9F, , , 1 ,,,/ , cfv%. , .d 5 “ cus , 


<9u£ 


dr 

dg" dvq ,, d/4' 


(p,q)£S 


We now use (8.6) to equate the first term in the expression above. In view of (8.7), the second 
term cancels out and we obtain, for the particular choice of F = F r , r : 1... M: 

Vr E [0, 77 ] Vr : 1 ... M 


X X X] ^(9,h)f^{g,h,{vP(g,h)}) 

0=1 —n (p,g)es 9 \(l,t)£S t 

8v P \ 

+ -^L( g ,h)-f™(g,h,{v p q (g,h)}) j =0. 

We proceed by evaluating the obtained formula at r = 0. Since dg/dr / 0 is an arbitrary vector 
in M n , and since the matrix |£^(x, {v|}) is invertible by (8.3), it follows that: 


OF 


dva . 


Va : 1... n V(p, q) E 5 


( 8 . 8 ) 


E If* 1 - Pl»/7(*. Id), {<( 1 . Pi})) 




+ {«&) - #*(*, Mh {«?(*, {«&)}) = 0, 


for all x E 77 and {u|} E V. 

3. Consider the following system of total differential equations: 


(8.9) 


Va : 1... n V(Z, t) 0 S 


dw 




f = /a t ( a; >{« ; t}(M)0S>{ l7 ?( X >{ U, t})}(P,9)£s)- 


<9Xr 


We now verify the Thomas condition for the system (8.9). It requires [2T] the following expressions 
to be symmetric in j. k : 1... n, for all (/, s) 0 5 on 77 x V: 




+ E 


{p,q&S) 


(l,t)£S 


XX™ r„-.D — ' ^ 
dwa 


By (8.8), the above expression equals: 

ei.s 


*^(x, {»'},{<})+ £ ^(x,{»'},{<})/J''(x, {»'},{<}) 

J (M)0S ° Wt 

O f*’ S 

+ E pEp.iwD.KH/f’fx, {»'},{<}), 

(P,5)GS 9 
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where as usual Vq = Vg(x,w l f ). We see that the required symmetry follows exactly by (8.6). 

Consequently, the problem (8.9) with the initial condition w\(xq) = {wq)\ has a unique solution 
onU - 


4. Let now {R'g}(p, g )eS be dehned on U through the formula: Wq(x) = Vq(x, We 

will prove that: 


( 8 . 10 ) 


Va : 1... n V(p, q) 6 S 


dwq 

dx Q 


fa q (*> { W t}(l,t)?S, {uJq}(p,q)£S 


Together with (8.9), this will establish the result claimed in the theorem. 


Differentiating (8.4) where we set v = w, and applying (8.9) yields: 


Mr : 1... M Va : 1... n Vx E U 
d F r 
dwl 


dF r , . .. dF r , . . .. dF r , . dwq , 

(x,w{x)) + 2^ -^ T (x,w(x))f t i(x,w(x))+ 2^ oTp( x >^( x ’))w—(®) = 0. 


dx 


(IMS 


(p,q)eS 


In view of (8.6), applied to F = F r , r = 1...M and v = w, the desired e qua lity in (8.10) 


follows directly by the invertibility of the matrix 
sufficiently small neighborhood O of (xq,wq). 


dF r 

dw% 


, X , w 


(P.9)GS 


, which by (8.3) is valid in a 


Remark 8.2. Another similar approach [lj, would be as follows. In case the Thomas condition 


(8.1) is not satisfied, one could relax the initial condition (4.7) and instead add the set of equa¬ 


tions (4.6a) to those in (8.1), before proceeding to derive the collections of equations in (8.2) by 


successive differentiation and substitution, as described above. Assume then that (8.3) holds in a 
large domain rather than a given point (xo,wo). The advantage of this method is that one does 
not have to be concerned about satisfying the initial condition potentially limiting the choices of 


the functions F r in (8.3), and hence one has a larger set of equations among (8.2 ) to c hoose from. 
Through this method, existence of an n 2 — M parameter family of solutions to (4.6c) follows 


as 


one has the freedom to set up the initial data in step 3. The disadvantage is that a larger set of 


equations in (8.5) must be satisfied for the sufficient condition to hold true. 


Remark 8.3. While Theorem 8.1 does provide exact algebraic conditions of integrability for (1.3), 


these conditions are by no means easy to verify for specific problems. With a view towards a more 
practical, if only approximate, solution procedure, the following alternative may be considered. 
In the spirit of Section [2j that sets up an infinite-dimensional optimization problem whose zero- 
cost solutions correspond to exact solutions of (1.3) (and whose non-zero cost solutions may 
be considered as approximate solutions of ( |1.3[ )), we now construct a family of finite-dimensional 
optimization problems parametrized by r G 11 (uncoupled from one x to another), whose pointwise 
zero-cost solutions (for each x) form the ingredients of an exact solution of (1.3), as outlined in 
Section [6j As before, we use the Einstein summation convention for all lowercase Latin indices 
from 1 to n. 
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Let w E M nxn , i £ !] C 1" and define: 

Aij(w,x ) = Gij(x) - w t i Gst(x)w s j 

= w s m f%(x) - ^G sm (x) ( wfdjGmpix ) + w‘-d i G mp (x) - w p j w^d t G pq (x)W^ l (w 

QJ3— f)B s OB— 0B S 

Cti k (vJ,x) = - -^£(w,x) + '-(w,x)B r pk (w,x ) - 




dwl 


dwl 


The symmetry of Bfj in the lower indices and C = 0 represent the integrability conditions of 


(4.6b) and (4.6c), while A = 0 is equivalent to (4.6a). Define: 


( 8 . 11 ) 


£{w;x) = \Aij(w,x)\ 2 + K 2 ^ I Cf jk (w,x)\‘ 


i,j =1 


i,j,k,s =1 


with Ki > 0 ,i = 1, 2 as nondimensional constants, where we assume that A. C in (8.11) have been 
appropriately nondimensionalized. 

We now seek minimizers w{x) of £{■ \x) for each x E fb Each set {u>(x); x E D} such that 
£(w(x)\x) = 0 on fl, defines a function x i-A w(x), ie!l, Each such function can be checked to 
see if it satisfies (4.6c), and any that passes this test allows for the construction of a solution to 


(4.6b) (because of the symmetry in i,j ) that in turn is a solution of (1.3) in view of A = 0. 

If the system (4.6b) and ( 4.6c| ) were completely integrable (i.e. satis fying the Thomas condition) 
then there would be a (?r 2 + n)-parameter family of solutions to (1.3). While the situation here is 
likely to be more constrained (with non-existence expected in the generic case), it is not possible 
to rule out, a-priori , the existence of solution families with fewer parameters. Thus, it is natural 
to expect nonuniqueness in looking for minimizers of £(■ ; x). This idea has the merit of lending 
itself naturally to a computational algorithm. 
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